The chir ally-odd twist-3 distribution function e a (x) 
in the chiral quark-soliton model 
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The chirally-odd twist-3 nucleon distribution e a (x) is studied in the large-JV c limit in the frame- 
work of the chiral quark-soliton model at a low normalization point of about 0.6 GeV. The remark- 
able result is that in the model e a (x) contains a 5- function-type singularity at x = 0. The regular 
part of e a (x) is found to be sizeable at the low scale of the model and in qualitative agreement with 
bag model calculations. 
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INTRODUCTION 



In deeply inelastic scattering (DIS) processes the nucleon structure up to twist-3 is described by six parton 
distribution functions, the twist- 2 fx(x), <7i(x), hf(x) and the twist-3 e a (x), g^(x), h%{x). Among these 
{SJ ' functions the least considered one is probably the twist-3 chirally odd distribution function e a (x) [1, 2], 

| which contrasts the fact that it is related to several interesting phenomena. E.g., a known [3-7] but rarely 

emphasized [8] fact is that e a (x) contains a (5-function-type singularity at x = 0, as follows from the QCD 
equations of motion. The existence of a <5(x)-contribution also was concluded from pcrturbative calculations 
[9]. The first Mellin moment of e a {x) is due to the ^(x)-contribution only. Another interesting phenomenon 
is connected to the first Mellin moment of the flavour-singlet (e u + e d )(x) which is proportional to the 
pion- nucleon sigma-term cr xN . The latter gives rise to the so-called "sigma-term puzzle". The large value 
extracted from pion-nucleon scattering data, cr xN « (50 — 70) MeV [10, 11], implies that about 20% of the 
nucleon mass M N is due to the strange quark - an unexpectedly large number from the point of view of the 
OZI-rule. 

The reason why e a {x) has received only little attention so far is related to its chiral odd nature, which 
means that e a (x) can enter an observable only in connection with another chirally odd distribution or 
fragmentation function, and therefore is difficult to access in experiments. Only recently it was shown that 
e a (x) can be accessed by means of the "Collins effect" [12], i.e. the left-right asymmetry in the fragmentation 
of a transversely polarized quark into a pion. This effect is described by the chirally and T-odd fragmentation 
function H^~ a (z), which is "twist-2" in the sense that its contribution to an observable is not power suppressed 
[12, 13]. First experimental indications to H± were reported in [14]. Assuming factorization it was shown 
that the Collins effect gives rise to a specific azimuthal (with respect to the axis defined by the exchanged hard 
virtual photon) distribution of pions produced in DIS of longitudinally polarized electrons off an unpolarized 
proton target. The observable single (beam) spin asymmetry is proportional to e\ e a (x) H^ a {z) [15] 
(e a = ±|, ±i are the quark electric charges). The process was studied in the HERMES experiment and the 
effect found consistent with zero within error bars [16]. 1 However, in the CLAS experiment, in a different 
kinematics, a sizeable asymmetry was observed [17, 18]. If the interpretation applies, that the CLAS data 
[17, 18] are due to the Collins effect, then e a (x) is definitely not small. Using estimates of H^(z) from 
HERMES data [19] it was shown in [20], that e a (x) could be about half the magnitude of the unpolarized 
twist-2 distribution /f (x) at Q 2 ~ 1.5 GeV 2 in the region 0.15 < x < 0.4 covered in the CLAS experiment. 

The indication that e a (x) could be large in the valence- a; region is not surprising, if one considers results 
from the bag model [2, 21], the only model where e a (x) has been studied so far. In this note e a (x) will be 
studied in the chiral quark-soliton model (xQSM). A subtle question is whether models with no gluon degrees 
of freedom (bag model, xQSM) can describe twist-3 distribution functions. The answer given in [1, 2] is yes, 
because e a (x), g^{x) and h a L {x) are special cases of more general quark-gluon-quark correlation functions; 
special inasmuch they do not contain explicit gluon fields. However, implicitly gluons do contribute and it 
is important and instructive to carefully interpret the results. E.g., in the bag model e a (x) is due to the bag 



1 The prominent result of the HERMES experiment [16] is the observation of sizeable azimuthal asymmetries in pion production 
from DIS of unpolarized electrons off longitudinally polarized protons, which contain information on H^ a (z) and the chirally 
odd distribution functions h^(x) and h^(x) [15]. 
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boundary [2]. This can be understood considering that the bag boundary (in a most intuitive way) models 
confinement, and thus "mimics" gluons. 

The xQSM was derived from the instanton model of the QCD vacuum. An important small parameter 
in this derivation is the "instanton packing fraction" which characterizes the diluteness of the instanton 
medium [22, 23]. Gluon degrees of freedom appear only at next-to- leading order of this parameter [24]. 
In leading order of the instanton packing fraction the xQSM quark degrees of freedom can be identified 
with the QCD quark degrees of freedom. This allows to consistently describe twist-2 quark and antiquark 
distribution functions f\(x), g"(x) and h1{x) at a low scale around 600 MeV [25]. The numerical results 
[25-29] agree within (10 — 30)% with parameterizations for the "known" distribution functions performed 
at low scales [30]. 

The twist-3 distribution functions g^(x) and h1(x) were studied in the instanton vacuum model in 
Refs. [31, 32]. The remarkable conclusion was that the pure twist-3 interaction dependent parts g^{x) 
and h>1(x) in the Wandzura-Wilczek(-like) decompositions of g^(x) and h1(x) are strongly suppressed by 
powers of the instanton packing fraction. As g^{x) and h1(x) do not contain explicit gluon fields it is possible 
to evaluate g?(x) and h1(x) directly in the xQSM. This was done in Refs. [33, 34], and it was observed that 
the pure twist-3 parts g^(x) and h^(x) are indeed small. Thus in these cases the xQSM respects the results 
found directly in the instanton vacuum model - i.e. the theory from which it was derived. This experience 
encourages to also tackle the study of e a (x) in the xQSM. However, one has to keep in mind that the results 
and their interpretation presented here should also be reexamined in the instanton vacuum model. This is 
out of the scope of this note and left for future studies. 

The xQSM describes the nucleon as a chiral soliton of the pion field in the limit of a large number of 
colours N c . This note focuses on the flavour singlet distributions (e u + e d ){x) and (e u + e d )(x) - the leading 
flavour combinations in the large--/V c limit. The consistency of the approach is checked, by demonstrating 
that the model expressions for (e u + e d )(x) satisfy QCD sum rules. It is shown that the model expressions are 
(quadratically and logarithmically) UV-divergent, and a consistent regularization is defined. Remarkably, it 
is found that the xQSM-expression for (e u + e d ){x) contains a (J(x)-contribution. The UV-behaviour and the 
<5(x)-contribution make an exact numerical evaluation of (e u + e d )(x) and (e u + e d )(x) particularly involved. 
Therefore (e u + e d )(x) is evaluated using an approximation, the "interpolation formula" of Ref. [25]. 

This note is organized as follows. In Sec. II the twist-3 distribution function e a (x) and some of its 
properties are discussed. Sec. Ill contains a brief introduction into the xQSM. In Sec. IV the flavour-singlet 
distribution (e u + e d )(x) is discussed and evaluated in the model using the interpolation formula. In Sec. V 
e a (x) is discussed in the non-relativistic limit. Sec. VI contains a summary and conclusions. 



II. THE DISTRIBUTION FUNCTION e a (x) 

The chirally odd twist-3 distribution functions e q (x) for quarks of flavour q and e q (x) for antiquarks of 
flavour q are defined as [1, 2] 

e q (x) = J_ y g e iA* {N \ ^(o) [0, \n] i> q (\n) \N) , e*» = e"(-x) , (1) 

where [0, An] denotes the gauge- link. The scale dependence is not indicated for brevity. The light- like 
vectors n a in Eq. (1) and p a are defined such that n a p a — 1 and the nucleon momentum is given by 
P£ = n a + p a M 2 /2. The matrix element in Eq. (1) is averaged over nucleon spin, i.e. (N\ . . . \N) = 
±Es 3 (N,S s \...\N,S 3 ). 

The renormalization scale evolution of e a (x) was studied in Refs. [3-5], see also Refs. [6, 7] for reviews. It 
evolves according to an evolution pattern typical for twist-3 quantities. It is not sufficient to know the n th 
Mellin-moment M n [e a ](Qo) = / dx x n ~ 1 e a (x, Q 2 ) at an initial scale Qq, in order to compute M n [e a ](Q 2 ) 
for Q 2 > Qq. Instead, the knowledge of all moments Mk[e a ](Qo) with fc < n is required. In the limit of a 
large number of colours N c the evolution of e a (x) simplifies to a DGLAP-type evolution - as it does for the 
other two proton twist-3 distributions h a L (x) and (the flavour non-singlet) g^(x). 

The QCD-equations of motion allow to decompose e g (x) in a gauge- invariant way as [3-5], see also [6, 7] 
and [8], 

e<1{x) = 6{X) 2^ W^ )^ )!^ + e q tw3 (x) + ^ ^ S n>1 . (2) 

The 5(x)-contribution has no partonic interpretation. Some authors cancel it out by multiplying e q (x) by x, 
while others prefer to consider from the beginning an alternative definition of e q (x) with an explicit factor of 
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x on the RHS of Eq. (1). The existence of a (^(^-contribution also was concluded in Ref. [9], where e q {x) was 
constructed explicitly for a one- loop dressed massive quark. The contribution e^ w3 (x) in Eq. (2) is a quark- 
gluon-quark correlation function, i.e. the actual "pure" twist-3 ("interaction dependent") contribution to 
e a (x), and has a partonic interpretation as an interference between scattering from a coherent quark-gluon 
pair and from a single quark [1, 2]. Its first two moments vanish. The third ("mass"-) term in Eq. (2) 
vanishes in the chiral limit. The "Kronecker symbol" S n> i accompanying it has the following meaning. The 
first moment of the mass-term vanishes. For x ^ the expression in Eq. (2) is, however, correct and can be 
used literally to take higher moments n > 1. 

The first and the second moment of e g (x) satisfy the sum rules [1, 2] 



dxeq{x) = 2j^< N \M°)M0)\N) , (3) 



dxxe^x) = ^N q . (4) 

In Eq. (4) N q denotes the number of the respective valence quarks (for proton N u — 2 and = 1). The 
sum rules (3, 4) follow immediately from the decomposition in Eq. (2) and the above mentioned properties of 
e tw3( x ) an d the mass-term. In particular, the sum rule (3) is saturated by the (5(x)-contribution in Eq. (2). 

The flavour singlet distribution function (e u + e d )(x) is related to the scalar isoscalar nucleon form- factor 
a(t) defined as 

a(t)u N (P')u N (P)=m(N(P')\ ^ u (0)^(0) + ^(0)^(0)) \N(P)) , t = (P' - P) 2 , (5) 

where un denotes the nucleon spinor (normalized as unUn = 2M N ) and m — \ (m u + md). In (5) a term 
proportional to (m u — md)('tpuipu — ^>dflpd) is neglected. At the point t = the form-factor a(t) is referred to 
as the pion-nucleon sigma-term ct tN and related to (e u + e d )(x) by means of the sum rule (3) as 



cr xN = cr(0) = m j dx(e u + e d )(x) . (6) 

The form factor a{t) describes the elastic scattering off the nucleon due to the exchange of a spin-zero particle 
and is not measured yet. However, low-energy theorems allow to deduce its value at the Cheng-Dashcn point 
t = 2mJ from pion-nucleon scattering data. One finds 

CT f2m 2 )-<f( 64±8 ) MeV Rcf ' [10] (7) 
a ^ m *> ~ \ (79 ± 7) MeV Rcf. [11]. (l) 

The difference cr(2mj) — er(0) was obtained from a dispersion relation analysis [35] and chiral perturbation 
theory calculations [36] with the consistent result of 14 MeV. This yields 

a„ N w (50 - 70) MeV . (8) 

With m w (5 — 8) MeV one obtains a large number for the first moment of (e u + e d )(x) 
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Ax {e u + e d ){x) w 10 . (9) 



One should keep in mind, however, that the large number in Eq. (9) is not due to a large "valence" structure 
in (e u + e d )(x), but solely due to the i5(x)-contribution. 

III. THE CHIRAL QUARK SOLITON MODEL (xQSM) 

The xQSM is based on the effective chiral relativistic quantum field theory given by the partition function 
[22, 37, 38] 




Z cii ■= VtpVtpVU exp i d^x ip(i@ -MU 1& - m)ip ) , (10) 



U = exp(iT a TT a ) , U^ 5 = ^(U + Ui) + \{U- ^ f )7 5 ■ (11) 
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In Eq. (10) M is the (generally momentum dependent) dynamical quark mass, which is due to spontaneous 
breakdown of chiral symmetry. U denotes the SU{2) chiral pion field. The current quark mass m in Eq. (10) 
explicitly breaks the chiral symmetry and can be set to zero in many applications. For certain quantities, 
however, it is convenient or even necessary to consider finite m. The effective theory (10) contains the Wess- 
Zumino term and the four-derivative Gasser-Leutwyler terms with correct coefficients. It has been derived 
from the instanton model of the QCD vacuum [23, 37], and is valid at low energies below a scale set by the 
inverse of the average instanton size 



600 MeV . 



(12) 



In practical calculations it is convenient to take the momentum dependent quark mass constant, i.e. M(p) — > 
M(0) = 350 MeV. In this case is to be understood as the cutoff, at which quark momenta have to be cut 
off within some appropriate regularization scheme. It is important to remark that (Mp av ) 2 is proportional 
to the parametrically small instanton packing fraction 



(Mp av ) 2 cx 



«1, 



(13) 



with i? av denoting the average distance between instantons. The smallness of this quantity has been used 
in the derivation of the effective theory (10) from the instanton vacuum model [23, 37]. 

The xQSM is an application of the effective theory (10) to the description of baryons [37]. The large- 
N c limit allows to solve the path integral over pion field configurations in Eq. (10) in the saddle-point 
approximation. In the leading order of the large- N c limit the pion field is static, and one can determine the 
spectrum of the one-particle Hamiltonian of the effective theory (10) 



H\n) = E n \n) , H = -i-y°^ k d k + 7 Mt/ 75 + 7 °m . 



(14) 



The spectrum consists of an upper and a lower Dirac continuum, distorted by the pion field as compared to 
continua of the free Dirac-Hamiltonian 



H \n ) = E no |n ) , H = -i^^dk + 7 °M + 7 °m , 



(15) 



and of a discrete bound state level of energy £q ov , if the pion field is strong enough. By occupying the 
discrete level and the states of the lower continuum each by iV c quarks in an anti-symmetric colour state, 
one obtains a state with unity baryon number. The soliton energy E so \ is a functional of the pion field 



E soi [U] = N c (E lcv + ( E n - Eno)) 

E n <0 



(16) 
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-E S oi[E^] is logarithmically divergent and has to be regularized appropriately, as indicated in (16). Minimiza- 
tion of E so i[U] determines the self-consistent soliton field U c . This procedure is performed for symmetry 
reasons in the so-called hedgehog ansatz 



7r a (x) = e a r P{r) , C/(x) = cos P(r) +ir a e* sin P{r) , 



(17) 



with r — |x| and e" = x a /r, in which the variational problem reduces to the determination of the self- 
consistent soliton profile P c {r). The nucleon mass M N is given by E so i[U c ]. The momentum and the spin and 
isospin quantum numbers of the baryon are described by considering zero modes of the soliton. Corrections 
in the l/./V c -expansion can be included by considering time dependent pion field fluctuations around the 
solitonic solution. A good and for many purposes sufficient approximation to the self-consistent profile P c (r) 
is given by the analytical "arctan-profile" 



( R 2 

P(r,m K ) — — 2arctan( 



■so/ 



(1 + m 7r r)e' 



Rsoi = M~ 



(18) 



where R so i is the soliton size. In the chiral limit m — > in (10) the self-consistent profile P c {r) cx 1/r 2 
for r — > co, and then — in (18). If m ^ in (10) the self-consistent profile exhibits a Yukawa-tail 
P c (r) cx e~ m7,r jr for r — > oo and is the physical pion mass connected to the current quark mass m in 
(10) by the Gell-Man-Oakes-Renner relation (see below Eq. (58)). The analytical profile (18) simulates this 
and the Gell-Mann-Oakes-Renner relation holds approximately. 
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The xQSM allows to evaluate in a parameter-free way nucleon matrix elements of QCD quark bilinear 
operators as (schematically) 

(N\^( Zl )T^(z 2 )\N) = c r 2M N A c ]T /d 3 X $ n ( Zl - X)r$„(z 2 - X) e iE ^-^ + ... (19) 

n occ 

= -c r 2M N A c /d 3 X $„( Zl - X)r$„(z 2 - X) e^"(^-^ ) + . . . . ( 2 0) 

n non 

In Eqs. (19, 20) T is some Dirac- and flavour-matrix, c r a constant depending on T and the spin and flavour 
quantum numbers of the nucleon state \N) = | S"3 , T3 ) , and <&„(x) = (x|n) are the coordinate space wave- 
functions of the single quark states \n) in (14). The sum in Eq. (19) goes over occupied levels n (i.e. n with 
E n < E\ cv ), and vacuum subtraction is implied for E n < E\ cv analogue to Eq. (16). The sum in Eq. (20) goes 
over non-occupied levels n (i.e. n with E n > £q ov ), and vacuum subtraction is implied for all E n > E\ cv . 2 
The dots in Eqs. (19, 20) denote terms subleading in the l/A c -expansion, which will not be needed in this 
work. Depending on the Dirac- and flavour-structure the expressions (19, 20) can possibly be UV-divergent 
and need to be regularized. If in QCD the quantity on the LHS of Eq. (19) is normalization scale dependent, 
the model results refer to a scale roughly set by in Eq. (12). 

In the way sketched in (19, 20) static nucleon properties (form-factors, axial properties, etc., see [39] for 
a review), twist-2 [25-29] and twist-3 [33, 34] quark and antiquark distribution functions, and off-forward 
distribution functions [40] have been studied in the xQSM. As far as those quantities are known, the model 
results agree within (10 — 30)% with experiment or phenomenology. It is important to note the theoretical 
consistency of the approach, in particular the quark and antiquark distribution functions in the model satisfy 
all general QCD requirements (sum rules, positivity, inequalities, etc.). 



IV. e a (x) IN THE xQSM 

A. Expressions and consistency 

a. Model expressions. The model expressions for the flavour combinations (e" ± e d )(x) "follow" from 
the expressions for the unpolarized twist-2 distributions (/" ± ff){x) derived in Ref. [25] by "replacing" 
the relevant Dirac-structure n a j a in the definition of fi(x) by 1/M N . This can be checked by an explicit 
calculation which closely follows the derivation given in [25] and can therefore be skipped here. This "anal- 
ogy" between e a (x) and ff(x) is due to the fact that both are "spin average" distributions, and the relevant 
Dirac- and flavour-structures exhibit the same properties under the hedgehog symmetry transformations. 
As a consequence, the flavour combinations (e u ± e d )(x) have the same large- N c behaviour as (/" ± ff)(x) 
[25], namely 

(e u + e d )(x) = N 2 c d(N c x) 

(e u -e d )(x) = N c d(N c x) , (21) 

where the functions d(y) are stable in the limit N c — > 00 for fixed arguments y = N c x, and different for the 
different flavour combinations. Though derived in the xQSM, the relations in (21) are of general character, 
considering that the xQSM is a particular realization of the large- picture of the nucleon [41] . The relations 
(21) are already the end of the story of "analogies" between e a (x) and fi(x) in a relativistic model. In the 
non-relativistic limit, however, e a (x) and fi(x) become equal, see Sec. V below. 

In this work only the leading order in the large-A c limit will be considered. At this order the model 
expressions read 

(e u + e d )(x) = A C M N ^ (n\-f°S(xM N - p 3 - E n )\n) (22) 

n occ 

= - A C M N (n\l°S(xM N -f- E n )\n) (23) 

n non 

and (e u — e d )(x) = as anticipated in (21). 



2 The possibility of computing model expressions in the two ways - (19) or (20) - has a deep connection to the analyticity 
and locality properties of the model [25]. In practice it provides a powerful check of numerical results. 
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b. Sum rule for the first moment. The first moment of the model expression in Eq. (22) reads 

1 

/ da; (e« + e d )(x) = N C Y] (n| 7 °|n) = ^ . (24) 
J ^ m 

_^ nocc 

When integrating over x in (24) one can substitute x — > y = xM N and extend the y-integration range 
[— M N ,M N ] to the whole y-axis in the large iV c -limit. The final step in (24) follows by recognizing in the 
intermediate step in (24) the model expression for the scalar isoscalar form- factor a(t) 

a(t) = mN c /d 3 x J (V~t |x|) Y, $ r;(x)7°*n(x) , (25) 

n occ 

at t = 0. (The Bessel-function Jq(z) = —* 1 for z — > 0.) The pion-nucleon sigma-term cr, N = <t(0) was 
studied in the xQSM in [42] and the form- factor <r(t) in [43]. 

The model expression for a„ N can also be derived in an alternative way using the Feynman-Hellmann 
theorem (this method was used in [42]) 

dM N (m) 

a„ N = m - B — . (26) 
Rewriting the expression for the nucleon mass M N = M N (m) in Eq. (16) as 

M N (m) = N c Y E n = N c Y (n\H\n) = M N (0) + m iV c £ (n| 7 °|n) , (27) 

nocc nocc nocc 

where vacuum subtraction is implied, and inserting M N (m) (27) into (26) one recovers the model expressions 
for a xN in Eqs. (24, 25). This proof is formally correct but one should be careful about regularization. A 
comment on that will be made at the end of Section IV C. 

c. Sum rule for the second moment. The second moment of (e u + e d ){x) in Eq. (22) is 



l 

dx x (e» + e d ){x) = ^]T (n| 7 V + E n )\n) - £ £7„ (n| 7 °|n) , (28) 



nocc nocc 



where (n| 7 °p 3 |n) drops out due to hedgehog symmetry. In QCD the sum rule (4) follows from using equations 
of motion. In the model the analogon is to use E n \n) = H\n). One obtains 

l 

/N 
dx x (e u + e d ) (x) = — - (m + 0M) , (29) 
M N 

-l 

^EH^W. (30) 

n occ 

where the relation E n (n| 7 °|n) = \ (n\{H, 7 °}|n) = m + M(n\ \{U + W)\n) was used. Eq. (29) then follows 

from X)nocc( n l n ) = Sp[O(£q cv + — H) — 0(— Ho)] — B where the vacuum subtraction is considered explicitly 
and B = 1 denotes the baryon number [25]. (Sp is the functional trace which can be saturated by respectively 

S P[-"] =E„allH ■■■!«> OT E„ all( n o|--.|«o).) 

For (3 = the QCD sum rule (4) would hold "literally" in the xQSM. However, in the model the equations 
of motions are modified compared to QCD and one cannot expect (3 = in (29, 30). Instead, the modified 
equations of motion in the xQSM suggest to interpret [3M (in the chiral limit) as the effective mass of model 
quarks bound in the soliton field. (One cannot expect (3 = 1 either, which would imply an effective mass 
M. It is jargon to refer to M as mass, strictly speaking M is a dimensionful coupling of the fermion fields 
to the chiral background field U .) 



B. Calculation of (e u + e d )(x) 



d. Interpolation formula. The approximation referred to as interpolation formula [25] consists in exactly 
evaluating the contribution from the discrete level to (e u +e d )(x) in Eq. (22), and in estimating the continuum 
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contribution as follows. One rewrites the continuum contribution in terms of the Feynman propagator in the 
static background soliton field U and expands it in powers of gradients of the [/-field, keeping the leading 
term(s) only 3 . The interpolation formula yields exact results in three limiting cases: (i) low momenta, 
| Vf | <c M, (ii) large momenta, \VU\ 3> M, (iii) any momenta but small pion field, | log E7| <C 1. One can 
expect that it yields useful estimates also in the general case. Indeed, it has been observed that estimates 
based on the interpolation formula agree with results from exact (and numerically much more involved) 
calculations within 10% [25-27]. 

e. Discrete level contribution. The Hamiltonian H (14) commutes with the parity operator tt and the 
grand-spin operator K, defined as the sum of the total quark angular momentum and isospin operator. The 
discrete level occurs in the = + sector of the Hamiltonian H (14). In the notation of Ref. [25] the 
discrete level contribution reads 

oo 

(e u + e d )(x) lcv = N C M N (lev| 7 °<5(xM N - E lcv - p 3 )|lcv) = N C M N J ^ (h(k) 2 - j(k) 2 ^j , (31) 

|xM N -B lev | 

where h(k) and j(k) are the radial parts of respectively the upper and lower component of the discrete level 
wave function in momentum-space, $i ev (p) = (p|lev), see [25] for details. 

/. Dirac continuum contribution. The two equivalent expressions, Eqs. (22, 23), allow to compute the 
contribution of the continuum states to (e u + e d ){x) in two different ways 

(e" + e d )(a;) cont = N C M N ^ (n\^°S{xM N - f - E n )\n) (32) 

E n <0 

= - N C M N {n\l°5{xM K -f- E n )\n) . (33) 



E n >0 



The expressions (32) and (33) can be rewritten by means of the Feynman propagator in the static background 
pion field as (see [25]) 



(e^ + e^cont = f>" + e d )(x)^ t , 

n=0 

(e u + e d )(x)£l - ImJV c M N J-^ S(p° + p 3 - X M K ) tr (p\(/>- UV > 



1 



iMWU-* |p) 



pi + V 2 - M 2 + iO V pi + V 2 - M 2 + iO 



(34) 



The vacuum subtraction is implied in (34) and means that the same expression but with U — ► 1 has to be 
subtracted, i.e. it is relevant only for the case n = 0. The subscript rog reminds that the expression might 
be UV-divergent and has to be regularized appropriately. Closing the p°-integration contour to the upper 
half of the complex p°-plane yields (32). Closing it to the lower half-plane yields (33). 

g. Gradient expansion: Zeroth order. Performing the traces over 7-matrices in the zeroth order contri- 
bution to (e u + e d )(x) cont in (34) yields 

(e u + e d )(x)^ nt = B sol A(x) (35) 

with 

B sol = ^ / d 3 xtr F [ ^— — - 1 I (36) 




A(x) = N C M N 8M Im / 7 ^ J ^ \^ ■ , n K/ . (37) 

rog 



It should be noted that this is not a strict expansion in gradients of the [/-field. The dimensionless parameter characterizing 
this expansion is l/(MR so i). Since the soliton solution is given for MR ao i = 1, see Eq. (18), such a strict expansion is 
not defined [25]. In the following higher orders in the gradient expansion will be considered merely in order to study the 
C/V-behaviour of the model expressions. 
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The coefficient B so i is real and contains the information on the soliton structure. The "—1" under the 
flavour-trace in (36) is due to vacuum subtraction. The function A(x) in (37) is well-defined within some 
appropriate regularization scheme to be figured out in the following. A(x) is an even function of x, provided 
the regularization is consistent with the substitutions p° — > — p° and p 3 — ► — p 3 in (37). Keeping x ^ and 
integrating over p 3 and p° (with the above-mentioned prescription to close the contour) one obtains 



^(•^)non/occ 



6(±x) 2N C M 



(2n) 2 



(38) 



rcg 



with the "+" sign referring to (33) and the "— " sign referring to (32). Thus A{x) is quadratically divergent, 
and depends on whether one computes it by means of (32) or (33). The non-equivalence of the two ways to 
compute a quantity in the model, (32) and (33), at the level of unregularizcd model expressions is a known 
phenomenon [25, 29]. The equivalence of (32) and (33) - and more generally of (19) and (20) - is a basic 
property of the model (see footnote 2). Therefore it is necessary to restore the equivalence of (32) and (33) 
in the expression (38) by means of a suitably chosen regularization. A regularization - which does this - is 
a Pauli-Villars subtraction of the type 



A(x)=A(x,M)-^-A(x,M 1 ) , 



(39) 



where Mi > M is the Pauli-Villars mass. The Pauli-Villars subtraction is the privileged method to regu- 
larize divergent distribution functions in the xQSM. This regularization preserves all general properties of 
distribution functions (QCD sum rules, positivity, etc.) [25], and where necessary it restores the equivalence 
of (32) and (33) in the final regularized model expressions [25, 29]. In the regularization (39) - which is 
sufficient at this stage - both formulae (32) and (33) yield the same result for A(x), namely 



A(x) = for x ^ 0. 



(40) 



Instead of studying the function A{x) at the point x = it is more convenient to consider moments of 
A(x) defined as .M„[A] = J dx x n ^ 1 A{x). Since A(x) is even in x, one needs to consider only odd moments 
2fc + l (k = 0,1,2,...) 



M[A]W+V = -^-Im 



K 

EC 

3=0 



d i p (p ) 2 '^) 2 * -2 ' 



(2tt) 4 p\ 



p 2 — M 2 + i0 



(41) 



rcg 



Performing a Wick rotation in Eq. (41), which is well defined for all moments under the regularization (39), 
one obtains 



M[A] 



(2fc+l) _ 



8iV c M 
M 2k 



E(S)(-i)' 



3=0 



d 4 p E (p° E ) 2j ( P l) 2k - 2j 



(42) 



Using 4-D spherical coordinates p£ = q (cos tf> , sin ^ sin 9 cos < 
k = q 2 one has 



(2tt) 4 p 2 + M 2 
, sin tp sin 9 sin <p , sin xp cos 9) and substituting 



M[A] 



(2fe+l) = _ 



8N c Ma k fdKn k+1 



f dfi 

2k J K- 



2(2ir) 4 M 2k J k + M 2 



rcg 



2tt 



a k = J2 (if) Jd<j>Jd9 JdiP sin6» sin 2 V (cosip) 2j (sin %j) cos9) 2k - 2j = 2n 2 5 k0 . (43) 

J =0 

I.e. only k = contributes which means that only the first moment of A(x) is non-zero 



M[A] 



(i) 



N C M 
2 ir 2 



dn k 



K+M 2 



M[A] in) = Ofor n > 2. 



(44) 



rcg 



The regularization prescription (39) removes the leading quadratic divergence in the integral in Eq. (44), 
but leaves a logarithmic one unregularized. However, e.g., a twofold Pauli-Villars subtraction 



A(x) = A(x, M) - aiA(x, Mi) - a 2 A(x, M 2 ) 



(45) 
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with 

M Mo — M 2 M Ml-M 2 

"1 = TT" 779 , "2 = -— — | — ~ , M 2 > M\ > M (46) 

Mi Mf - M 2 M 2 M| - Ml ' 

is sufficient to make the first moment of M[A]^^ finite. In the limit Mi — > oo one recovers the previous 
regularization (39). It is important to note that the modification (45, 46) of the previous regularization (39) 
still preserves the property (40). Thus, (e u + e d )(x)fj nt (x) satisfies 

(e» + e d ) (*)<!(*) - for x ? 0, | dx x n ~\e u + e d ) = C S nl , (47) 

i.e. (e M + e d )(x)fj nt (x) is proportional to a (5-function at x = 0. What remains to be done is to compute the 
coefficient C of the <5- function in Eq. (47). 

The two Pauli-Villars subtractions in Eq. (45) introduce two parameters, Mi and M 2 in Eq. (46), which 
have to be fixed. E.g. one could first fix Mi by regularizing the logarithmically UV-divergent model 
expression for the pion decay constant f„ 



U C J (27T 



d 4 p E M 2 



)4 (p| + M 2 f 



(48) 

rog 



such that it gives the experimental value f„ = 93MeV. Then one could fix M 2 , e.g., by means of the quark 
vacuum condensate which is given in the effective theory (10) by the quadratically divergent expression [39] 



(vac\(ip u ip u + ipd-ipd^vac) = -8N C 



d 4 p E M 



(2tt)4 pi + M 2 



(49) 

rog 



Two subtractions analogue to (45) are required to regularize (49). In this way the free parameters Mi and 
M 2 are fixed. 4 It is important to stress that the parameters Mi and M 2 are fixed in the vacuum- and in 
the meson-sector of the effective theory (10). In this sense the xQSM, i.e. the baryon sector of the effective 
theory (10), produces parameter-free results. 

It is interesting to observe that the coefficient C can directly be expressed in terms of the quark condensate 
(49), such that 

C = B so i (vac|(^ u V« + VWd)|vac) (50) 

with B so i defined in (36). In Eq. (50) any details on the regularization have disappeared. However, the 
explicit demonstration of the existence of a regularization prescription, which preserves the properties in 
Eq. (47), was a crucial step in the derivation of (50). 

h. Gradient expansion: Higher orders. Taking the trace over 7- and flavour-matrices in the expression 
for the first order term in the expansion (34), (e" + e d )(x)^J nt , one obtains 



(e" + e^X^Ocont — ImiV c M N (—i8M) J ^ ^ /(p| V fc cosP(f, m v )\p) 



0, (51) 



since (plV* cosP(f,m 7r )|p) = J d 3 xV fe cos P(r, m ff ) = 0. A lengthy calculation yields for the second order 
contribution in (34) the result 

(e" + e<*)(4L = Im^g^ /| j da J d P fd^tr F A 



A = U\x-ave 3 )[\7 k U(x- f3ve 3 )][\7 k u\x)} + (U ^ f/ f ) , (52) 



4 For f% in (48) a single subtraction, f% = fl{M) - {M/M{) 2 }1 (Mi) with Mi = 556 MeV is required. For the quark 
condensate (49) one needs two subtractions analogue to (45, 46) with M2 = (9.1 ± 5.7) GeV in order to reproduce the 
phenomenological value [(280 ± 30) MeV] 3 [44]. The first Pauli-Villars mass Mi is of the order of magnitude of the "natural 
cutoff" pav 1 ai 600 MeV of the effective theory (10). The much larger second Pauli-Villars mass M2 - in some sense merely 
introduced as a "technical device" to remove the "residual (logarithmic) divergence" left after the subtraction of the leading 
(quadratic) divergence - has no physical meaning. In non-renormalizable theories (such as the xQSM) the cutoff (here 
Mi ~ pay 1 ) has a physical meaning and shows up in final expressions. In renormalizable theories the result for the quadratically 
divergent integral (44) would be proportional to M 2 , see e.g. [45]. 
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FIG. 1: a. The oversimplified partonic interpretation of e q (x). b. The symbolic diagrammatic interpretation of 
the 8(x) contribution to e a (x) as suggested by Eq. (54). The quark-lines "carry the fraction x = of the nucleon 
momentum", see text. 

which does not depend on how contours are closed, i.e. the two ways, (32) and (33), yield the same result. 
In (52) terms are neglected which contain three or more derivatives acting on the [/-fields. In a strict 
expansion in the number n of [/-field gradients these terms have to be considered in higher orders n > 3. 
The result (52) is real and UV-finite. Still, it has to be regularized according to the prescription (45, 46). 
In an exact evaluation of the continuum contribution it would be, of course, not possible to pick up the 
divergent contribution from the zeroth order in the gradient expansion (34) and regularize only that. Since 
(e u + e d )(x)l?J nt oc M, the application of the regularization (45, 46) to (52) yields 

(e" + e rf )(x)fl = 0. (53) 

All higher orders n > 3 in the gradient expansion (34) are UV-finite. 

i. Intermediate summary. In [29] it was shown that - if it occurs - the non-equivalence of (32) and 
(33) at the level of unregularized model expressions only shows up in the lowest UV-divergent order(s) of 
the expansion (34). Thus, the results of this section show that the regularization (45, 46) (i) consistently 
regularizes the continuum contribution to (e u + e d )(x), and (ii) ensures the equivalence of (32) and (33). 

The final regularized model expression for the continuum contribution consists of a <5-function at x = 
with a coefficient proportional to the quark condensate and the factor B so i in Eq. (36) which encodes the 
information on the nucleon (i.e. soliton) structure 

(e" + e d )(a;) con t = C 5(x) , C = B sol (vac|(^„V« + VWd)|vac) . (54) 

The existence of the 5- function is a feature of the model (with the Pauli-Villars regularization method) . 
The fact that the continuum contribution consists of no regular part but the (5-function only has to be 
considered as a peculiarity of the approximation (interpolation formula) used. 

j. How to interpret a S-function? A (5(x)-contribution to a distribution function has no partonic inter- 
pretation. The observation in (54), however, suggests an intuitively appealing "interpretation". 

Oversimplifyingly e a (x)dx can be interpreted as taking out of the nucleon in the infinite momentum 
frame, e.g., a left-handed good light-cone quark component which carries between x and x + dec of the 
nucleon momentum, and then reinserting a right-handed bad light-cone quark component with the same 
momentum back into the nucleon, see Fig. la. (The good and bad quark light-cone degrees of freedom are 
strictly speaking defined in the light-cone quantization. Good means independent degrees of freedom, the 
bad quark degrees of freedom are composites of good quark and gluon degrees of freedom. Considering this 
one obtains the correct partonic interpretation of e a {x) as a quark-gluon-quark correlation function [1, 2].) 

Does it make sense to pick up hereby a quark (or antiquark) which carries the fraction x — of the 
nucleon momentum, i.e. which is at rest with respect to the fast moving nucleon? Eq. (54) suggests that 
such a quark (or antiquark) is picked up from the vacuum, which to a certain extent is present also inside the 
nucleon. It should be stressed that one does not deal with a disconnected diagram. The factor B so i shows 
that the nucleon line and the vacuum blob are connected by the exchange of a resonance with the quantum 
numbers of the sigma meson, see the symbolic diagram in Fig. lb. 

It would be interesting to see whether such an "interpretation" could be confirmed by observations analogue 
to (54) in other models. 
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FIG. 2: a. The regular part of the chirally odd twist-3 flavour-singlet quark and antiquark distribution functions, 
(e u + e d )(x) and (e u + e d )(x), at a low normalization point of the model of about 600 MeV vs. x. At x = there is 
a 5(a;)-contribution. b. The chiral quark-soliton model results for (e u + e d )(x) (computed here) and (/" + f d )(x) 
(from [27]) at a low scale of about 600 MeV vs. x. c. The same as Fig. 2b but for antiquarks. 



C. Discussion of the results for (e™ + e d )(x) 

The final result for (e u + e d )(x) from the interpolation formula is the contribution of the discrete level 
(31), (already the total result for x ^ 0) and the continuum contribution (54) consisting of a <5(a;)-function. 
It should be noted that there is no freedom to also regularize the UV-finite discrete level contribution in the 
Pauli-Villars regularization method. This contribution must not be regularized fore otherwise the variational 
problem of minimizing the soliton energy in (16) has no solution, i.e. no soliton exists [27]. 

Fig. 2a shows the final results for (e u + e d )(x) and (e u + e d ) (x) (no effort is made to indicate the (5-function 
at x = 0). It is instructive to compare (e u + e d ){x) to (/" + f d )(x) in the model. Both are of the same 
order in the large iV c -limit and become equal in the non-relativistic limit (see Sec. V). For the quarks one 
observes that (/" + f d )(x) is about 2-3 times larger than (e u + e d )(x), while the corresponding antiquark 
distributions are of a similar magnitude, see Figs. 2b and 2c. 

In order to compute the coefficient C of the <5(a;)-function in Eq. (54), one has to evaluate B so i (36). For 
the physical situation with m w = 140 MeV one obtains 

B sol = ^ d 3 xtr F (^±^ - l) = / d3x (cosP(r,m T ) - lj = -17.2 R 3 sol . (55) 

In the chiral limit — > the integral in the expression for B so i in (36) can be evaluated in an elementary 
way yielding B so i = —2^/2 1: 2 R? sol = —27.9 Rl ol - I.e. the coefficient C is by about 40% increased in the chiral 
limit. This demonstrates the importance of considering this quantity in the physical situation with finite 
m n . Numerically one obtains for R so i — M~ x = (350 MeV) -1 and = 140 MeV the result 

C= (9.1 ±2.8), (56) 

where the error is due to the uncertainty of the vacuum condensate (—280 ± 30) 3 MeV 3 [44]. For the first 
moment of (e u + e d )(x) one thus obtains 

l l l 

J Ax (e u + e d )(x) = Jdx {e u + e d )(x) lev + jdx {e u + e d )(x) cont = 1.6 + (9.1 ± 2.8) = (10.7 ± 2.8) (57) 

-l -l -l 

in good agreement with (9). In order to obtain cr xN from Eq. (57) it is convenient to use the Gell-Mann- 
Oakes-Renner relation 

m lfl = - m (vac|(^„V« + VWd)|vac) , (58) 
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which holds in the effective theory (10) [39] and allows to eliminate the uncertainty from the phenomenological 
value of the vacuum condensate in the continuum contribution to a wN 

K N )cont = mj &x{e u + e d )(x) cont =mC= -m%fi B sol = 67.8 MeV. (59) 

In order to obtain the contribution of the discrete level to a^ N one can use (58) to obtain a consistent value 
for the current quark mass m = (8.2 ± 2.5) MeV. This yields (cr„ N ) lev = (13.2 ± 4.1) MeV and the total 
result is (j 7tN — (81 i 4) MeV. There is no point in keeping track of the error in this case, since it is smaller 
than the accuracy of the interpolation formula (which was found to be about ±10% whenever it was checked 
quantitatively). Thus one obtains 

a„ N = 81 MeV . (60) 

The result (59) is about 30% larger than former exact results from the xQSM which, however, have been 
calculated with a different (proper-time) regularization [42, 43]. Considering that cr xN is quadratically di- 
vergent and thus rather strongly sensitive to regularization, the result in (60) is in good agreement with 
the results of Rcfs. [42, 43]. Worthwhile mentioning is that all model numbers - from Eq. (60) and from 
Refs. [42, 43] - are consistent with the phenomenological value for cr xN in Eq. (8) within (10 — 30)%. 

In QCD - as mentioned in Sec. II - the first moment of e a {x) is due to the <5(x)-function only. In the xQSM 
the (5(x)-function provides the dominant (more than 80%) but not the only contribution to the first Mcllin 
moment of (e" + e d )(x). The second moment of (e u + e d )(x) receives no contribution from the continuum 
and is due to the discrete level contribution only 

J dxx (e" + e d ) (x) = 0.20 . (61) 

This result would imply that the quarks bound in the soliton field have an effective mass f3M ~ 78 MeV (in 
the chiral limit), see the discussion below (29, 30). 

A comment is in order on an exact numerical evaluation of (e u + e d )(x). The distribution functions 
computed in the xQSM so far were all either UV-finite or at most logarithmically divergent. In the latter 
case always a single Pauli-Villars subtraction was sufficient. The regularization prescription (45, 46) precisely 
states how (e u + e d )(x) can practically be regularized. However, an exact numerical evaluation meets the 
problem to evaluate the model expressions for a Pauli-Villars mass M-i k, several GeV (see footnote 4) In 
the numerical calculation the spectrum of the Hamiltonian (14) is discretized and made finite (see e.g. [26] 
and references therein). For the latter step one considers only quark momenta below some large numerical 
cutoff A num chosen much larger than any other (physical or numerical) scale involved in the problem. So 
far A num w (8 — 10) GeV was sufficient, but this is the order of magnitude of the second Pauli-Villars mass 
M 2 . To compute (e u + e d )(x) one has to choose A num much larger than M 2 , which would result in an 
uneconomically large increase of computing time. 

Interestingly, the singular (5(x)-contribution would conceptually cause no problem for the numerical method 
of Ref. [26]. The descretized spectrum of the Hamiltonian (14) yields discontinuous (distribution) functions 
of x. In [26] it was proposed to smear the distribution functions q(x), i.e. to convolute them with a narrow 

Gaussian with an appropriately chosen width 7 as g(a;) smcar = ^ l/2 J dx' exp(— ^ x ~2 ^ )q( x ')- (The smearing 

can be removed by a deconvolution procedure.) This trick would turn the (5(x)-contribution in (e u + e d )(x) 
into a narrow-Gausssian with the well-defined width 7. In this way the coefficient of the 5(x)-contribution 
could be well determined from the numerical result. 

Finally, a comment is in order on (26) which relates the logarithmically divergent nucleon mass M N and 
the quadratically divergent cr„ K . M N requires a single Pauli-Villars subtraction, while a wN requires two 
subtractions. Thus the quantities on the left-hand and right-hand side in (26) are regularized differently 
In the Pauli-Villars regularization scheme the relation (26) has to be considered as formally correct modulo 
regularization effects. (Some other ambiguities in the Pauli-Villars regularization scheme were mentioned 
in [46].) In other regularization methods - such as the proper time regularization - there are no such 
ambiguities. 

D. Comparison to the bag model 

Studies of e a (x) were also performed in the framework of the bag model in Refs. [2, 21]. In Fig. 3 the 
xQSM result for the regular part of (e u + e d )(x) is compared to the result from the MIT bag model from 
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FIG. 3: The twist-3 distribution function (e u + e d )(x) from the xQSM (solid line, obtained here) and from the bag 
model (dashed line, from Ref. [2]) vs. x. The meaning of the curves at negative x is explained in Eq. (1). The xQSM 
result refers to a scale of about 0.6 GeV. The bag model result refers to a scale of about 0.4 GeV. 



Rcf. [2]. (For that the flavour- independent results for "e(x)" from Refs. [2] are multiplied by the factor 
N c = 3 in order to compare to (e u + e d )(x) obtained here.) The comparison of the position of the maxima 
of the curves for quark distributions from the two models indicates that the bag model results refer to a 
somehow lower scale than 0.6 GeV, the scale of the xQSM. (In [21] the value of 0.4 GeV was quoted.) Taking 
this into account one concludes that both models give qualitatively similar results for quark distributions. 

Concerning antiquarks the difference is more pronounced. However, the bag model description of antiquark 
distributions cannot be considered as reliable. A drawback of the bag model in this context is that it yields 
negative /f (x) in contradiction to the positivity requirement. 



V. e a (x) IN THE NON-RELATIVISTIC LIMIT 

In the limit of the soliton size R so i — ► the expressions of the xQSM go into the results of the non- 
relativistic ( "naive" ) quark model formulated for an arbitrary number of colours N c [47] . In this sense the 
limit R so i — ► corresponds to the non-relativistic limit in the xQSM. This was studied in detail in Ref. [48]. 

As R so i — ► the soliton profile (18) goes to zero, and the [/-field approaches unity. Correspondingly, 
the spectrum of the Hamiltonian (14) becomes more and more similar to that of the free Hamiltonian (15). 
Considering vacuum subtraction it is clear that the contribution of the continuum vanishes in this limit [48] , 
and all that remains is the contribution of the discrete level |lev). More precisely, as R so i — * the energy 
of the discrete level E\ cv — > M (so the nucleon mass, (16), formally M N — > N C E\ CV — > N C M), and the lower 
component of the Dirac-spinor of the discrete level wave- function goes to zero [48] . 

The limit R so i — -> means that U — ► 1 i.e. log U <C 1 which is the case (iii) in which the interpolation 
formula yields exact results (see Section IV B). The first feature - in this case the vanishing of the continuum 
contribution in (32, 33) - can be observed in the final result (e u + e d )(x) cont = C S(x) in Eq. (54). The factor 
B so i in the coefficient C vanishes with R so i — > as can be seen from its definition (36) or Eq. (55). Thus, in 
the non-relativistic limit the contribution of the 5-function at x = vanishes because the coefficient C goes 
to zero. 

To study the non-relativistic limit in the discrete level contribution it is convenient to use the first ex- 
pression in Eq. (31). Since only the upper component of the Dirac-spinor of the discrete level wave- function 
survives the limit [48], one can replace 7 by the unity matrix, i.e. 

lim (e u + e d )(x) — > N c M N (\ev\S(xM N -f- £W)|lev) . (62) 

Next consider that E\ ev — ► M and M N — » N C M while the momenta of the non-relativistic quarks \p l \ -C M 
such that the corresponding operator in the (5-function in (62) can be neglected. Using the normalization 
(lev I lev) = 1 one obtains 

lim (e u + e d ){x) — » N c s( x - . (63) 
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The first two moments of (63) read 

1 l 
lim [dx (e u + e d )(x) = N c , lim [dx x (e u + e d ){x) = 1 . (64) 

non-rel / non-rel / 

-l 



To see that the relations (64) are the correct non-relativistic results for the QCD sum rules (4) and (6) one has 
to consider that in the non-relativistic limit the current quark mass m — > M — M N /N C and a„ N — > M N . (The 
latter relation follows formally, e.g., from the Feynman-Hcllmann theorem (26) with m — M — M N /N C .) 
The twist- 2 unpolarized flavour-singlet distribution function (/" + ff)(x) is given in the xQSM by [25] 

(fx + /i )(*) = N C M N H(l + 7°7 3 ) S(xM N -p 3 - E n )\n) , (65) 

n occ 

i.e. the only difference to (e u + e d )(x) is the different Dirac-structure (1 + 7°7 3 ) instead of 7 . This difference 
becomes irrelevant in the non-relativistic limit (where only the upper component of the Dirac-spinor of the 
discrete level wave- function survives) such that in this limit (e u + e d )(x) and (/" + f d )(x) become equal. 
This argument holds also for separat flavours and allows to generalize 

lim e\x) = lim f q (x) = N q s(x - -M, (66) 

non-rel non-rel y J\ Q J 

where N q denotes the number of the respective valence quarks (i.e. for the proton N u = 2 and Nd = 1 for 
N c = 3 colours). The result (66) means that in the non-relativistic limit e q {x) is given by the mass term 
contribution in the decomposition (2) because 

-£ *,.(&>]- to/jw. (e7, 



i[M N x J M N 



n-rol 



In the intermediate step in (67) the S(x — 1/N C ) -function was used to replace x in the denominator by l/N c . 
One may worry that the large value cr xN = M N would mean a large strangeness content of the nucleon. 
However, as discussed in [8] the value <t„ n = M u correctly implies a vanishing strangeness contribution to 
the nucleon mass in the non-relativistic limit. 

Thus, though phenomenologically it is not satisfactory, the non-relativistic picture of the twist-3 distribu- 
tion function e q (x) is consistent. In this limit the singular and pure twist-3 contributions in the decomposition 
(2) vanish, and e q (x) is given by the mass-term. Thus the chirally odd nature of e q {x) arises from a "mass 
insertion" into a quark line. Moreover, e q (x) and j\ (x) become equal in this limit, and are trivial (5-functions 
concentrated at x = 1/N C which means that the nucleon momentum is distributed equally among the N c 
massive and non- interacting constituent quarks. 

The usefulness of results of the kind (66) is best illustrated by the popularity of the non-relativistic relation 
between the twist- 2 helicity gf(x) and transversity h\(x) distribution functions (here for N c = 3) 

lim h\{x)= lim g q (x)=P q s(x-\) , P u = \ , P d = - \ , (68) 

non-rel non-rel \ O / O O 

which yields for the axial charges g 1 ^ = | and g^ = 1. Though also these numbers are phenomenologically 
not fully satisfactory the assumption that h\(x) = g1{x) at some low scale is a popular guess to estimate 
effects of transversity distribution, see [49] for a review. 



VI. SUMMARY AND CONCLUSIONS 



A study of the flavour-singlet twist-3 distribution function (e u + e d ){x) in the xQSM was presented. 
It was shown that the model expressions are quadratically and logarithmically UV-divergent and can be 
regularized by the Pauli-Villars method. The model expressions for the quark and antiquark distribution 
functions (e u + e d ){x) and (e u + e d ){x) were evaluated using an approximation - the interpolation formula 
which in general well approximates exact model calculations. 

The remarkable result is that the xQSM-expression for (e u + e d )(x) contains a (5-function-type singularity 
at x — as expected from QCD [8]. This result is obtained here from a non-perturbative model calculation. 
Previously a <5(x)-contribution in e a (x) was observed in a perturbative calculation in Rcf. [9]. 
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In the xQSM the coefficient of the (5-function is proportional to the quark vacuum condensate. This 
is natural from the point of view that in QCD the singular contribution to e q (x) and the quark vacuum 
condensate are both the expectation values of the same local scalar operator, ip q (0)ip q (0), taken respectively 
in the nucleon and vacuum states. This observation allows to make a heuristic but physically appealing 
interpretation of the (5(x)-contribution. 

At x 7^ the xQSM yields results for (e u + e d ){x) similar to those obtained in the bag model at a 
comparably low scale [2, 21]. Both models suggest that e a (x) is sizable at low scales. The discription of 
(e u + e d )(x) in the xQSM is consistent in the sense that the sum rules for the first and the second moment are 
satisfied. However, in the xQSM the 5(x)-contribution provides the dominant but not the only contribution 
to the first moment of (e u + e d )(x) unlike in QCD, and in the case of the second moment it is necessary 
to interpret the result correspondingly by introducing the notion of an effective quark mass. It would be 
interesting to see whether the failure of the bag model to satisfy the sum rule for the second moment reported 
in [2] could also be reinterpreted in a similar spirit. 

In effective models, such as xQSM or bag model, equations of motions are altered compared to QCD and 
there is no gauge principle which would allow to cleanly decompose e a {x) into a <5(x)-contribution, a pure 
twist-3 part and a mass-term. Therefore one cannot expect that the QCD sum rules which are derived by 
means of the QCD equations of motion are literally satisfied in such models. Still within the models the 
results are consistent. 

The (large- N c ) non-relativistic limit of e q (x) was studied on the basis of the xQSM expressions. It was 
found that in this limit e q (x) — ff{x). The non-relativistic description of e q (x) was shown to be consistent. 

The results and interpretations presented here should be reexamined in the more general framework of 
the instanton model of the QCD vacuum, on which the xQSM is founded. This was out of the scope of the 
study presented here and will be reported elsewhere. 
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Note added in Proof: After this work has been completed the work [50] appeared, where the authors con- 
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[1] R. L. Jaffe and X. D. Ji, Phys. Rev. Lett. 67, 552 (1991). 
[2] R. L. Jaffe and X. D. Ji, Nucl. Phys. B 375, 527 (1992). 

[3] I. I. Balitsky, V. M. Braun, Y. Koike and K. Tanaka, Phys. Rev. Lett. 77, 3078 (1996). 
[4] A. V. Belitsky and D. Muller, Nucl. Phys. B 503, 279 (1997). 
[5] Y. Koike and N. Nishiyama, Phys. Rev. D 55, 3068 (1997). 

[6] A. V. Belitsky, in Proceedings of the "31st PNPI Winter School on Nuclear and Particle Physics", St. Petersburg, 

Russia, 24 Feb - 2 Mar 1997, ed. V. A. Gordeev, pp. 369-455 [arXiv:hep-ph/9703432]. 
[7] J. Kodaira and K. Tanaka, Prog. Theor. Phys. 101, 191 (1999). 
[8] A. V. Elremov and P. Schweitzer, arXiv:hep-ph/0212044. 
[9] M. Burkardt and Y. Koike, Nucl. Phys. B 632, 311 (2002). 
[10] R. Koch, Z. Phys. C 15, 161 (1982). 

[11] M. M. Pavan, I. I. Strakovsky, R. L. Workman and R. A. Arndt, ttN Newsletter 16, 110 (2002) [arXiv:hep- 

ph/0111066]. 
[12] J. C. Collins, Nucl. Phys. B 396, 161 (1993). 

X. Artru and J. C. Collins, Z. Phys. C 69, 277 (1996). 
[13] D. Boer and P. J. Mulders, Phys. Rev. D57, 5780 (1998). 

D. Boer, R. Jakob and P. J. Mulders, Phys. Lett. B424, 143 (1998). 

D. Boer and R. Tangerman, Phys. Lett. B381, 305 (1996). 
[14] A. V. Efremov, O. G. Smirnova and L. G. Tkachev, Nucl. Phys. Proc. Suppl. 74, 49 (1999); Nucl. Phys. Proc. 

Suppl. 79, 554 (1999). A. V. Efremov et al, Czech. J. Phys. 49, S75 (1999) [arXiv:hep-ph/9901216]. 
[15] P. J. Mulders and R. D. Tangerman, Nucl. Phys. B461, 197 (1996) [Erratum-ibid. B484, 538 (1996)]. 



16 



[16] A. Airapetian et al. [HERMES Collaboration], Phys. Rev. Lett. 84, 4047 (2000); H. Avakian [HERMES Collab- 
oration], Nucl. Phys. Proc. Suppl. 79, 523 (1999); A. Airapetian et al. [HERMES Collaboration], Phys. Rev. D 
64, 097101 (2001). 

[17] H. Avakian [CLAS Collaboration], talk at 9th International Conference on the Structure of Baryons (Baryons 

2002), Newport News, Virginia, 3-8 March 2002. 
[18] H. Avakian et al. [CLAS Collaboration], arXiv:hep-ex/0301005. 

[19] A. V. Efremov, K. Goeke and P. Schweitzer, Phys. Lett. B 522, 37 (2001) [Erratum- ibid. B 544, 389 (2002)]. 
[20] A. V. Efremov, K. Goeke and P. Schweitzer, Acta Phys. Polon. B 33, 3755 (2002); and arXiv:hep-ph/0208124. 
[21] A. I. Signal, Nucl. Phys. B 497, 415 (1997). 

[22] D. I. Diakonov and V. Y. Petrov, Nucl. Phys. B 272, 457 (1986). 
[23] D. I. Diakonov and V. Y. Petrov, Nucl. Phys. B 245, 259 (1984). 
[24] D. I. Diakonov, M. V. Polyakov and C. Weiss, Nucl. Phys. B 461, 539 (1996). 

[25] D. I. Diakonov, V. Petrov, P. Pobylitsa, M. V. Polyakov and C. Weiss, Nucl. Phys. B 480, 341 (1996). 

[26] D. I. Diakonov, V. Y. Petrov, P. V. Pobylitsa, M. V. Polyakov and C. Weiss, Phys. Rev. D 56, 4069 (1997). 

[27] C. Weiss and K. Goeke, arXiv:hep-ph/9712447. 

[28] P. V. Pobylitsa and M. V. Polyakov, Phys. Lett. B 389, 350 (1996); 

P. V. Pobylitsa et al., Phys. Rev. D 59, 034024 (1999); 

M. Wakamatsu and T. Kubota, Phys. Rev. D 60, 034020 (1999); 

P. Schweitzer et al, Phys. Rev. D 64, 034013 (2001). 
[29] K. Goeke, P. V. Pobylitsa, M. V. Polyakov, P. Schweitzer and D. Urbano, Acta Phys. Polon. B 32, 1201 (2001). 
[30] M. Gluck, E. Reya and A. Vogt, Z. Phys. C 67, 433 (1995). M. Gliick, E. Reya, M. Stratmann and W. Vogelsang, 

Phys. Rev. D 53, 4775 (1996). 
[31] J. Balla, M. V. Polyakov and C. Weiss, Nucl. Phys. B 510, 327 (1998). 
[32] B. Dressier and M. V. Polyakov, Phys. Rev. D61, 097501 (2000). 
[33] M. Wakamatsu, Phys. Lett. B 487, 118 (2000). 
[34] M. Wakamatsu, Phys. Lett. B 509, 59 (2001). 

[35] J. Gasser, H. Leutwyler and M. E. Sainio, Phys. Lett. B 253, 252 and 260 (1991). 
[36] T. Becher and H. Leutwyler, Eur. Phys. J. C 9, 643 (1999). 

[37] D. I. Diakonov, V. Y. Petrov and P. V. Pobylitsa, Nucl. Phys. B 306, 809 (1988). 

D. I. Diakonov and V. Y. Petrov, JETP Lett. 43, 75 (1986) [Pisma Zh. Eksp. Teor. Fiz. 43, 57 (1986)]. 
[38] D. I. Diakonov and M. I. Eides, JETP Lett. 38, 433 (1983) [Pisma Zh. Eksp. Teor. Fiz. 38, 358 (1983)]; A. Dhar, 

R. Shankar and S. R. Wadia, Phys. Rev. D 31, 3256 (1985). 
[39] C. V. Christov et al, Prog. Part. Nucl. Phys. 37, 91 (1996). 

[40] V. Y. Petrov, P. V. Pobylitsa, M. V. Polyakov, I. Bornig, K. Goeke and C. Weiss, Phys. Rev. D 57, 4325 (1998); 

M. Penttinen, M. V. Polyakov and K. Goeke, Phys. Rev. D 62, 014024 (2000); 

P. Schweitzer, S. Boffi and M. Radici, Phys. Rev. D 66, 114004 (2002). 
[41] E. Witten, Nucl. Phys. B 160, 57 (1979); E. Witten, Nucl. Phys. B 223, 433 (1983). 
[42] D. I. Diakonov, V. Y. Petrov and M. Praszalowicz, Nucl. Phys. B 323, 53 (1989). 
[43] H. C. Kim, A. Blotz, C. Schneider and K. Goeke, Nucl. Phys. A 596, 415 (1996). 
[44] J. Gasser and H. Leutwyler, Phys. Rept. 87, 77 (1982). 
[45] R. Delbourgo and M. D. Scadron, Mod. Phys. Lett. A 17, 209 (2002). 
[46] T. Kubota, M. Wakamatsu and T. Watabe, Phys. Rev. D 60, 014016 (1999). 
[47] G. Karl and J. E. Paton, Phys. Rev. D 30, 238 (1984). 
[48] M. Praszalowicz, A. Blotz and K. Goeke, Phys. Lett. B 354, 415 (1995). 
[49] V. Barone, A. Drago and P. G. Ratcliffe, Phys. Rept. 359, 1 (2002). 
[50] M. Wakamatsu and Y. Ohnishi, preprint OU-HEP-433, arXiv:hep-ph/0303007 (2003). 



